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ABSTRACT 


High  speed  oscillatory  flows,  like  high  speed  mean  flows,  are  capable  of  inducing 
time-averaged  heat  transfer  effects.  This  research  involves  the  analytical  solution  of  a  model 
problem  of  zero-mean  internal  oscillatory  flow,  which  arises  from  a  high-intensity  resonant 
standing  acoustic  wave  set  up  across  the  ends  of  two  parallel  plates.  The  compressible  form 
of  the  Navier-Stokes  equations  are  solved,  along  with  the  equations  of  continuity,  energy, 
and  state,  using  perturbation  solution  and  complex  variable  methods.  MAPLE,  a  symbolic 
mathematical  software  tool,  is  utilized  to  find  the  time-averaged  portion  of  the  temperature 
distribution  between  the  plates.  The  final  heat  transfer  results  are  presented  in  terms  of 
suitably  defined  recovery  factors.  The  analysis  is  performed  for  varying  gap  widths  between 
the  plates  using  air  as  the  host  fluid.  This  work  provides  the  fundamental  explanation  of  the 
phenomenon  responsible  for  the  thermoacoustic  refrigerating  effect  as  well  as  an  analytical 
basis  for  determining  the  optimum  gap  width  between  the  plates  of  the  stack  in  a 
thermoacoustic  refrigerator. 
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I.  INTRODUCTION 


This  research  involves  the  study  of  the  time-averaged  heat  transfer  effects  induced 
by  high  speed  internal  zero-mean  oscillatory  flows.  The  model  problem  being  treated  is 
for  an  ideal-gas  host  fluid  between  two  “infinite”,  parallel  plates.  The  zero-mean 
oscillatory  flow  is  produced  by  a  high-intensity  standing  axial  acoustic  wave  induced 
across  the  ends  of  the  plates.  The  plates  are  separated  by  a  distance,  lyo,  and  the  sound 
wave  has  a  wavenumber  of  k^lnlX  where  A,  is  the  wavelength.  The  relationship  between 
these  two  parameters,  yo  and  k,  is  assumed  to  be  such  that  yo^«l  throughout  the  analysis. 
In  other  words,  the  narrow  acoustic  waveguide  approximation,  which  is  discussed  further 
in  Chapter  HI,  holds  true  for  this  research.  Another  general  assumption  made  for  this 
analysis  is  that  the  host  fluid  between  the  plates  is  a  gas  which  is  treated  as  ideal. 

The  particular  focus  of  this  study  is  to  solve  analytically  the  Navier-Stokes 
equations  for  this  problem.  The  governing  equations  are  set  up  and  simplified  using  the 
particular  geometry  and  assumptions  pertinent  to  this  particular  case.  Then,  the  variables 
are  scaled  in  a  standard  fashion  to  reduce  the  equations  of  continuity,  momentum,  energy 
and  state  into  dimensionless  forms.  An  order  of  magnitude  analysis  is  performed  to 
further  simplify  the  equations.  Finally,  perturbation  methods  are  employed,  and  like- 
order  terms  are  grouped  to  form  new  equations.  The  equations  of  most  significant  order 
are  solved  analytically  using  complex  variable  methods  for  their  respective  unknowns. 
The  final  goal  is  to  determine  the  resulting  time-averaged  temperature  distribution  and 
study  its  affects  on  the  heat  transfer  behavior  at  the  walls.  MAPLE,  a  symbolic 
computational  software  package  (see  Blachman  and  Mossinghoff  (1994)),  is  utilized  to 
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perform  the  extensive  mathematical  manipulations  that  are  required  in  this  analysis.  The 
results  of  the  time-averaged  thermal  behavior  that  is  ultimately  of  interest  are  presented, 
by  convention,  in  terms  of  suitably  defined  recovery  factors. 

The  entire  analysis  is  performed  for  two  separate  cases,  which  differ  based  on  the 
width  of  separation  between  the  plates,  2yo,  in  comparison  with  a  length  scale  called  the 
Stokes’s  boundary  layer  thickness,  8^.  The  Stokes’s  boundary  layer  thickness,  also  called 
the  viscous  penetration  depth,  is  defined  in  Schlichting  (1979)  in  the  discussion  of 
Stokes’s  second  problem  concerning  the  flow  near  an  oscillating  plate  as  the  thickness  of 
the  layer  of  fluid  which  is  carried  by  the  plate.  If  the  quantity,  yo/S^,  is  greater  than  or 
equal  to  a  certain  amount  (which  is  determined  in  this  research  and  discussed  in  Chapter 
rV),  the  spacing  between  the  plates  is  considered  to  be  a  “wide  gap.’’  The  first  and 
simpler  case  is  analyzed  using  a  wide  gap  geometry  and  its  respective  assumptions.  In 
the  second  case,  an  arbitrary  gap  width  is  considered  which  lifts  some  simplifying 
assumptions  so  that  it  includes  “narrow  gaps.”  The  arbitrary  gap  case  is  useful  for 
analysis  of  the  time-averaged  transport  behavior  in  the  stack  region  of  a  thennoacoustic 
refrigerator. 
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II.  BACKGROUND 


Since  the  time-averaged  heat  transfer  results  of  this  work  are  presented  in  terms  of 
conventionally  defined  recovery  factors,  it  is  useful  to  review  the  concept  of  a  recovery 
factor.  In  order  to  do  so,  the  application  of  recovery  factors  to  two  well  studied  flows, 
namely  Couette  and  Poiseuille  flows,  will  be  presented  here  as  adapted  from  Burmeister 
(1983). 

Couette  flow,  as  illustrated  in  Figure  2.1,  occurs  when  fluid  exists  between  two 
infinite  parallel  plates. 


Figure  2. 1  Couette  Flow.  After  Burmeister  (1983). 


One  plate  (the  upper  plate  here)  is  moving  with  a  velocity  Vo,  and  the  other 
(lower)  plate  is  stationary.  The  plates  are  at  different  temperatures,  Jo  and  7,.  If  both 
plates  were  stationary,  the  temperature  distribution  between  them  would  simply  be  linear; 
however,  the  fluid  motion  caused  by  the  motion  of  the  upper  plate  affects  the  temperature 
distribution.  This  effect  comes  from  viscous  dissipation  and  is  especially  significant  in 
high  speed  flows.  Applying  the  following  assumptions  and  boundary  conditions  to  the 
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governing  equations  of  continuity,  momentum,  and  energy,  the  velocity  and  temperature 
distributions  are  determined. 

The  assumptions  being  made  are:  (1)  steady  conditions,  (2)  laminar  flow,  (3) 

constant  properties,  (4)  no  pressure  gradient  in  the  x  direction,  (5)  no  edge  effects  or  end 

effects,  (6)  Newtonian  fluid,  and  (7)  no  buoyancy  or  other  body  forces. 

Boundary  Conditions:  At  y=0:  u=0=w  (no  slip  at  wall) 

v=0  (no  penetration  at  wall) 

T=To  (fixed  wall  temperature) 

At  y=2yo:  u=Vo,  w=0  (no  slip  at  wall) 

v=0  (no  penetration  at  wall) 

T=Ti  (fixed  wall  temperature) 

The  resulting  velocity  distribution  is  linear  as  given  below: 


Vo  2y, 

The  resulting  temperature  distribution  is  given  by 


T-T, 

T,-To 


+  (£cPr)(^)(l- 
2yo 


(2.1) 


(2.2) 


V'^ 

where  Ec  =  - — — ^ - ,  which  represents  a  comparison  of  viscous  dissipation  effects  to 

2Cp(T^-To) 
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where  T„aii=To  from  before.  Now  it  is  observed  that  the  heat  transfer  is  driven  solely  by 
the  temperature  difference,  T„airTa^,  where  T^,  is  not  constant  but  varies  depending  on  the 
particular  geometry  and  flow  conditions  being  studied. 


This  result  incorporates,  into  a  simple  form,  the  effects  of  viscous  dissipation  on 
heat  transfer  which  become  more  and  more  significant  as  Vq  and  the  kinetic  energy  of  the 
fluid  increase.  It  also  allows  the  heat  transfer  coefficient  to  be  expressed  in  a  more 
consistent  and  convenient  manner  as 
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The  heat  transfer  coefficient,  as  a  result,  is  never  negative;  and  the  Nusselt  number, 
Nu  = - ^ ,  conveniently  has  the  value  of  one. 

K 

From  the  expression  for  the  adiabatic  wall  temperature  in  Equation  2.5,  it  is 

clearly  observed  that  some  of  the  fluid’s  kinetic  energy,  — ,  is  being  converted  into 

2 

thermal  energy.  In  a  sense,  some  of  this  kinetic  energy  is  being  “recovered”  as  thermal 
energy  at  the  lower  adiabatic  wall  which  then  results  in  a  rise  in  temperature  of  that  wall. 
The  ratio  of  the  thermal  energy  recovered  at  an  adiabatic  surface  to  the  kinetic  energy  of 
the  fluid  is  called  the  recovery  factor,  r.  For  this  Couette  flow,  the  recovery  factor  is 

-Ti)  _ 


The  physical  significance  of  the  recovery  factor,  as  a  measure  of  viscous 
dissipation  effects,  applies  equally  well  for  other  types  of  flow  and  wall  conditions,  even 
for  situations  with  non-adiabatic  surfaces.  For  any  situation,  “r”  contains  information 
about  how  much  of  the  fluid’s  kinetic  energy  is  converted  into  thermal  energy  at  the 
walls.  It  is  defined  and  calculated,  however,  in  a  similar  manner  to  that  of  the  Couette 
flow  example  -  the  ratio  of  the  recovered  thermal  energy  at  an  adiabatic  surface  to  the 
fluid’s  kinetic  energy.  The  expression  of  the  kinetic  energy  may  differ  depending  on  the 


velocity  chosen  to  base  it  on,  but  the  temperature  difference  used  for  calculation  is  always 
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Of  course,  for  other  types  of  flow,  the  resulting  expression  for  the  recovery 
factor  may  be  more  complicated  than  simply  a  linear  Prandtl  number  dependence. 

Another  type  of  flow  is  Poiseuille  flow  which  is  of  direct  relevance  to  the  current 
work.  In  Poiseuille  flow,  pictured  in  Figure  2.2,  both  surfaces  are  stationary.  The 
assumptions  applied  to  this  case  are  the  same  as  for  the  Couette  flow  case  with  the 
exception  that  there  is  a  uniform  pressure  gradient  in  the  x  direction  which  drives  a  uni¬ 
directional  mean  flow.  Using  the  following  boundary  conditions,  the  velocity  and 
temperature  distributions  are  determined. 


Figure  2.2  Poiseuille  Flow.  After  Burmeister  (1983). 


Boundary  Conditions: 


At  y=-yo: 

u=0=w 

(no  slip  at  wall) 

v=0 

(no  penetration  at  wall) 

T=To 

(fixed  wall  temperature) 

At  y=+yo: 

u=0=w 

(no  slip  at  wall) 

v=0 

(no  penetration  at  wall) 

7’=r, 

(fixed  wall  temperature) 

The  resulting  velocity  distribution  is  given  by 


2\i  dx 


yo 


(2.9) 
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The  definition  of  average  velocity, 


2>o 


(2.10) 


gives 


3^  dx^ 


(2.11) 


Also  the  maximum  velocity,  which  occurs  at  the  centerline  Cy=0)  is  given  by 


.  -yo  (  dp  3 

"■  2  ^  2“"'’ 


(2.12) 


The  velocity  distribution  from  Equation  2.9  may  then  be  expressed  in  terms  of  either  the 
average  or  the  maximum  velocity  as 
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«=«m(l-"V)  or  M  =  -M„,(l-^) 

2  yl 


(2.13) 


The  resulting  temperature  distribution  for  this  Poiseuille  flow  is  given  by 


^=7’o+^4^a+— )+^«^(i-4) 
2  Jo  3*c  yl 


(2.14) 


From  Fourier’s  law,  the  wall  heat  flux  at  the  lower  wall  is  given  by 


-riy  =  -^o)  = - ;; - - - 

dy  2yo 


(2.15) 


Using  the  concepts  established  in  the  Couette  flow  case,  the  adiabatic  wall  temperature 
for  Poiseuille  flow  is  deduced  to  be 


Taw 


6PrM 


(2.16) 
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This  can  be  also  easily  verified  by  solving  the  adiabatic  wall  case  and  finding  the 
temperature  of  the  adiabatic  wall. 

Again,  it  is  observed  that  the  driving  potential  for  heat  transfer  between  the  walls 
is  the  difference  in  temperature,  T^airTaw-  This  ensures  that  viscous  dissipation  effects  are 

]c 

accounted  for  and  that  the  heat  transfer  coefficient,  h  =  —^,  is  always  positive.  The 

2yo 


Nusselt  number,  Nu  =  ,  is  again  unity.  The  recovery  factor,  r,  is  then  given  by 


.  _  p 


cAT^-T,) 


=  12Pr 


u 


av 

2 


(2.17) 


Note  that  the  kinetic  energy  of  the  fluid  was  chosen  to  be  represented  in  terms  of  the 
average  velocity. 

These  concepts  of  an  adiabatic  wall  temperature  and  a  corresponding  recovery 
factor  apply  to  other  geometries  and  flow  conditions  and  are  typically  of  consequence  in 
high  speed  flows  where  velocity  gradients  and  viscous  dissipation  effects  are  significant. 
In  high  speed  mean  flows,  a  sufficient  amount  of  the  fluid’s  kinetic  energy  is  transformed 
into  thermal  energy,  causing  significant  heat  transfer  to  occur  at  the  walls  from  viscous 
dissipation.  Normally,  the  result  is  a  heating  of  the  walls. 

Recovery  factors  will  be  found  for  the  geometiy  and  flow  conditions  present  in 
this  research  which  are  similar  to  the  example  case  of  Poiseuille  flow  just  considered. 
The  physical  geometry  is  exactly  the  same  with  the  presence  of  two  parallel  plates.  The 
principal  difference  in  the  current  problem  is  in  the  flow  conditions.  In  this  case  a  driving 
(zero-mean)  oscillatory  pressure  gradient  is  set  up  across  the  plates  by  a  high-intensity 
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resonant  standing  acoustic  wave.  This  contrasts  to  the  uni-directional  mean  flow  which 
occurs  in  Poiseuille  flow.  Despite  being  a  zero-mean  oscillatory  flow,  the  cross-coupling 
between  the  harmonics  is  capable  of  inducing  a  time-averaged  flow  and  a  thermal  effect 
which  can  be  significant  in  high  speed  situations.  This  key  idea  introduces  a  number  of 
interesting  new  features  into  the  problem  which  are  explored  in  this  research.  The  posed 
problem  and  its  solution  will  be  addressed  further  in  Chapter  in,  and  the  corresponding 
recovery  factor  will  be  discussed  extensively  in  Chapter  IV. 
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III.  THEORY 


A.  THE  MODEL  PROBLEM  AND  GOVERNING  EQUATIONS 

The  model  problem  consists  of  high  speed  oscillatory  flow  between  two  “infinite” 
parallel  plates  and  is  depicted  in  Figure  3. 1 .  Both  plates  are  held  at  a  fixed  temperature, 
T„.  The  host  fluid,  air,  oscillates  between  the  plates  in  the  jc*  direction  with  a  velocity 
amplitude  of  oscillation,  Vo.  and  a  frequency,  o),  in  radians. 


Figure  3.1  Model  Problem:  Oscillatory  Flow  Between  Parallel  Plates. 


This  flow  is  induced  by  a  high  intensity  standing  acoustic  wave  set  up  across  the 
ends  of  the  plates.  The  acoustic  wave  number  is  h=2ii/\  where  A,  is  the  acoustic 
wavelength.  The  wavelengths  of  concern  in  this  research  are  on  the  order  of  1  meter. 
The  distance  between  the  plates  is  2yo,  and  a  Cartesian  coordinate  system  is  established 
with  the  origin  at  the  centerline.  Using  a  wavelength  of  the  above  mentioned  order,  1 
meter,  the  narrow  acoustic  waveguide  approximation,  yok«l,  requires  that  y©  be  much 
less  than  1  meter.  In  this  research,  the  gap  widths  (2yo)  studied  are  on  the  order  of 
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millimeters  or  tenths  of  millimeters.  Therefore,  the  narrow  acoustic  waveguide 
approximation  certainly  holds  throughout  this  work. 


To  determine  the  temperature  distribution  between  the  plates,  which  is  the 
fundamental  objective  of  this  study,  the  compressible  form  of  the  equations  of  continuity, 
motion,  and  energy  are  solved  analytically.  These  equations  in  Cartesian  coordinates  are 
(asterisks  (*)  denote  a  dimensional  quantity): 

Continuity:  ^  +  _(p  +^(p  * v*) +^(p  * w*)  =  0  (3.1) 

at*  dx*  ay*  az* 

Motion: 


X  *  component: 

9m*  *3w*  *  *  3|p* 

p  *(- —  +  u*- —  +  V*- —  +  H'*-^— )  =  p  *  2.,* — — 
dt*  9x  *  dy  *  dz  dx  * 

9  .^^u*  2  du*  9v*  9w*,, 

,  9  r.,^9M*  9v*  9  r  ,9>v*  9m*^, 


(3.2a) 


y  *  component: 

.*.dv*  9v* 


•  +  H'*— )  =  p*g,*-^ 


^  r  2  9m*  9v*  9w*^, 

9y*  ^  ^9y*  3 

9  9v*  9w*  9  r  /9m*  9v*.- 


(3.2b) 


z  *  component: 

p  ♦  - 1_  ^  - {_  Y  ♦  - \-w* - )  =  p  * 

dt*  dx*  dy*  dz  dz  * 

9  //,9w*  2  du*  9v*  9w*^_ 

9  9w*  9m*  9  r  /^v*  9w*.- 


(3.2c) 
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Energy: 

dr*  dr*  dr*  dr*  d  dr* 

P*C,(^  +  M*^+V*-^+W*^)  =  {— (it,— ) 
df*  dx*  dy*  dz  dx*  dx* 

^  dn  *  ^  dp  * 

+  V  *^f- +  w*-^)  +  jlO)  * 
dy  *  dz 


(3.3a) 


where:  4..  =  2[(|iil)= +(|::I)= +(^)=]  +  [^+i^l» 
dx  *  dy  *  dz  *  dx  *  dy  * 

^  dw*  ^  dv*  2  ^  r^M*  ^  dw*  2  2  dn*  dv*  dw*  2 

dy  *  dz  *  dz  *  dx  *  3  dx  *  dy  *  dz  * 


(3.3b) 


Note  that  in  the  energy  equation,  p  denotes  the  coefficient  of  thermal  expansion,  and  O* 
is  called  the  viscous  dissipation  function. 

Certain  basic  assumptions  being  applied  to  this  problem  are:  (1)  narrow  acoustic 
waveguide  approximation,  yo^«l,  (2)  no  spanwise  (z*)  dependence  and  no  edge  effects, 

=  0 ,  (3)  no  gravitational  field  effects,  g;,*=gy*=g^*=0,  (4)  no  imposed  heat 


generation,  q'"  =  0 ,  and  (5)  fluid  is  a  gas  which  may  be  treated  as  being  ideal  with 


constant  properties  (dynamic  viscosity,  p,  specific  heat  at  constant  pressure,  c^,  and 
thermal  conductivity,  The  following  boundary  conditions  are  also  imposed: 


At  y*  =  ±yo :  m*  =  0 

(no  slip  at  wall) 

V*  =  0 

(no  penetration  at  wall) 

II 

* 

(fixed  wall  temperature,  called  mean  temperature) 

du  * 

Aty*  =  0:  ”  =0 

dy* 

(symmetry) 

0 

II 

* 

> 

(“zero-crossover”) 

13 


The  governing  equations  are  transformed  into  dimensionless  equations  by  scaling 
the  dimensional  variables  using  the  standard  scales  for  the  flow  fluctuation  quantities  in 


an  acoustic  field.  The  scaling  relations  are: 


t  =  cot*,  X  = 


(kUn) 


u  *  >’  “  V' 

kx*,  y  =  — ,  M  =  — ,  v  = - 


(^3^0^) 


P-(P*-PJ  ,-P*-P'n  ^  = 

(P.M)  ’  ^  (p„V)’ 


where  the  Mach  number  is  Af  =  — ,  c  is  the  speed  of  sound,  p„  is  the  mean  mass  density. 


p„  is  the  mean  pressure,  and  the  reference  temperature,  AT„f,  is  as  yet  undefined. 

After  implementing  the  scaling  factors,  applying  the  assumptions,  and  utilizing 
algebraic  manipulations  in  order  to  simplify  the  expressions,  the  terms  with  order  of 
magnitude  of  1  are  neglected.  The  resulting  dimensionless  equations  are: 

Continuity:  m|-(pm)  +  ^+ Af:^(pv)+|^  =  0  (3.6) 

at  dx  dx  dy  dy 

X  -  Momentum: 

(pM  +  l)(— +  Mm— +Mv— )  =  -^  +  (— )—  ^  ^ 

at  dx  dy  dx  y^  dy 

where  the  boundary  layer  thickness  is  6  =  ^  =  ,  and  the  kinematic  viscosity  is 
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Using  order  of  magnitude  analysis  and  ylk^«  I,  the  pressure  gradient  term  is 
the  only  dominant  term  in  the  y-momentum  equation  (Equation  3.8)  which,  therefore, 

reduces  to  ^®0.  This  allows  the  boundary  layer  approximation  to  be  made  that 

dy 

pressure  is  dependent  only  on  x  and  t  (p  =  p(x,t)) . 

The  remaining  governing  equations  are  solved  utilizing  a  regular  perturbation 
analysis.  The  perturbation  method  assumes  a  solution  of  the  form: 

/  =  /i +e/2+e^/.+—  (3.10) 
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where  e  is  a  small  quantity.  In  this  problem,  the  appropriate  small  quantity,  the 
perturbation  parameter,  is  the  Mach  number,  M.  Therefore,  solutions  are  chosen  in  the 
forms; 


M  =  a,  +  Mmj  +  M  + ... 

=  T7‘l’o +<l>i  +  ^<l>2  +  +  ••• 

M 

p  =  p,  +  Mp2  +  M^p,  +  ... 
v  =  V,  +  A/Vj  +  +  ... 

p  =  p,  •  Mp2  +  MV3+  -- 


(3.11) 


Note  that  the  <|)  solution  includes  a  term  of  order  M  '  (written  0(M  ')),  since  a  term  of  that 
order  appears  in  the  energy  equation  (Equation  3.9). 

These  solution  forms  in  Equation  3.11  are  substituted  back  into  Equations  3.6, 
3.7,  3.8  and  3.9  which  are  subsequently  expanded  algebraically  and  ordered  according  to 
the  powers  of  M.  The  like  order  terms  are  grouped,  and  separate  sets  of  partial 
differential  equations  are  extracted  for  each  order.  The  most  significant  solutions  are 
those  with  orders  M  '  and  W  (or  0(1)),  with  the  remaining  higher  order  effects  not  being 
considered  in  this  study. 


Overall,  the  dimensionless  equations  to  be  solved  are 


Continuity: 


dp,  3m,  3v, 

^  +  ^  +  ^  =  0 
at  ax  ay 


(3.12) 


X  -Momentum: 


3m,  ^  I  (  J 
3t  dx  2r\l  3y^ 


(3.13) 
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where  r\^=-^  is  a  gap  width  parameter  based  on  the  viscous  penetration  depth, 

^  VJ 


Energy: 


O(M'): 


3/  2n5Pr  3y’  '  T,  ’ 


dt 


(3.14) 


Energy:  0(1): 


dt 


+  P 


a(|)o  ^  „  9<|) 


'  dt 


+  u, 


•  + V 


9(t) 


0  _ 


dx  ‘  dy  2r\l  Pr  dy 


1  3  ^(b ,  - 

- ^  +  2pr„Ec 


(3.15) 


The  variables,  u,  <|),  p,  v,  and  p  are  real  functions  of  time  as  well  as  spatial 
variables.  In  general,  a  real  variable,  Z,  can  be  expressed  as  the  real  part  of  a  complex 
function  as  follows: 


Z(x,y,0  =  Re[z(x,y)e‘“]  =  l[ze’"'  (3.16) 

where  z{x,y)  is  the  complex  amplitude,  z  **is  the  complex  conjugate  of  z(x,y),  and 
“n”  represents  the  harmonic.  Applying  this  concept  to  the  unknown  variables  in 
Equations  3.12, 3.13,  3.14,  and  3.15,  it  follows  that: 
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M,  (jc,  y,t)  =  Re[M,  (x,  y)e“  ]  =  +u^*  *e'‘'  ] 

<t)oU>>’.0  =  =  +4*0  **«""] 

4) ,  (x,  >, 0  =  <t>  ,0 (^,  y)  +  Re[^,2  (x,  y)e^‘'  ]  =  ([> ,o ( >)  +  :|- [^,2^'"  +  ^,2  *  *e"^"  ] 

1 

PiU,y,0  =  Re[p,(j:,yy']  =  -[p,e"  +p,  **e-"] 

Vi  (x, y,0  =  Re[ V,  (x, y)e“  ]  =  +  v,  *  *e'“  ] 

p, {x,t)  =  Re[p, (x)e“ ]  =  ^[p.e''  +  P,  *  *e~“ ] 

The  variables,  Ui,  ([Jo,  pi,  Vj,  and  pi  exhibit  leading  order  behavior,  whereas  <|>i,  Ui, 
P2,  and  all  others  contribute  to  higher  order  behavior.  These  two  categories  of  behavior 
are  analyzed  in  turn  in  the  following  sections.  The  behavior  of  the  pressure  distribution, 
however,  is  particularly  significant  to  this  work;  therefore,  the  analysis  of  p,  is  treated  in  a 
separate  section  following  the  discussion  of  the  other  leading  order  terms. 

C.  LEADING  ORDER  BEHAVIOR 


The  leading  order  equations  to  be  solved  are  as  follows: 

Continuity:  0(1):  =  0 

ot  dx  dy 


X  -Momentum:  0(1):  ^^  =  --^  +  (— ^)^ 

dt  dx  2r]ydy^ 

Energy:  O(M'):  %=._!„  +2|}r.Ec(l 


dt  2t1o  Pr  dy^ 


dt 


(3.18) 

(3.19) 

(3.20) 
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The  leading  order  equations  to  be  solved  are  now  Equations  3.22,  3.23,  and  3.24. 
Using  the  imposed  boundary  conditions  along  with  the  proper  equation  of  state, 
M, ,  ([) 0 ,  p , ,  V, ,  and  p,  are  found. 

The  equation  of  state  gives  density  as  a  function  of  pressure  and  temperature, 
p=p(p,7).  A  differential  element,  Ap*,  is  then  given  as 
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Now, 

(|^)r*  =  Ap*  =  p*-p„=  pp„V,c  (3.26) 

Also, 

Ar*  =  r*-r„  =0Ar„^  and  Ap*  =  p*-p„  =  PP„M  (3.27) 


From  the  definition  of  the  coefficient  of  thermal  expansion,  B  =  — (^^)  , ,  it  follows 

p*  dT*  ’’ 


(^)p-  =  -PP*  =  -Ppm(l  +  Mp)  (3.28) 

Utilizing  all  of  these  relations.  Equation  3.25  becomes 

PP.^  =  ;^PP«K,c-  P„p„(l  +  A/p)<|)Ar„^  (3.29) 

Substituting  in  the  perturbation  solution  forms  from  Equation  3.11,  the  following 
equation  results  for  pi: 

P.  ='lP.-P».7;(Y-l)<t>o  (3.30) 

Using  the  complex  notation  established  in  Equation  3.17,  it  results  that: 

P>='1P.-P„7;(Y-1)^o  (3.31) 

For  ease  of  following  this  discussion,  the  final  leading  order  equations  to  be  solved  and 
their  corresponding  boundary  conditions  are  summarized  below: 

^  ..  •  ^  diu, 

Continuity:  -^  =  -,p  ^  (3,32) 

oy  dx 
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(3.33) 

(3.34) 

(3.35) 


The  density,  p,,  comes  directly  from  the  solution  of  the  equation  of  state 
(Equation  3.35).  That  leaves  three  equations  and  four  unknowns  which,  at  first  glance, 
are  impossible  to  solve.  However,  a  fourth  equation  is  obtained  when  applying  one  of  the 
boundary  conditions  to  the  continuity  equation  (Equation  3.32).  When  applying  the 
condition  that  v,  =0  at  y=0,  an  equation  for  p^(x)  results.  Then,  using  the  appropriate 
boundary  conditions  for  p, ,  the  function,  p,(jc),  is  determined  explicitly.  The  expression 
for  p,(x)  is  then  utilized  in  the  remaining  equations  to  solve  explicitly  for  the  remaining 
unknowns,  m,  ,  <|)o,  p , ,  and  v,  . 

Using  standard  methods  of  solving  partial  differential  equations  (such  as 
separation  of  variables).  Equations  3.32,  3.33,  3.34,  and  3.35  were  solved.  Some 
simplifying  substitutions  were  also  made.  Since  the  reference  temperature  is  defined  to 
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be  ,  the  Eckert  number,  as  defined  in  Equation  3.9,  is  simply  unity.  Also,  for 

an  ideal  gas,  the  quantity,  is  approximately  unity.  Thus,  the  leading  order  solutions 
are  given  by  the  following  expressions: 


dx  cosh((l  +  i)riQ) 


(3.36) 


-r  .  .  .  .^1  cosh((l  +  i)T|oyVPr)- 

(t)o (X,  y)  =  2p,  (x)[l - —-7-.-'^'^  1 

cosh((l  +  j)rioyVPr) 


(3.37) 


cosh((l  +  j)TioVPr) 


(3.38) 


v,  (x, y)  =  lp^  (x)[(l  -  y)  +  (Y  - 1)( 


tanh((l  +  t)TioVPr) 

(l  +  0i1oV^ 


sinh((l  +  OqoyV]^)  _i) 

(l  +  i)T|oV^cosh((l+i)T|oV^)  dx^  ^ 

+  (tanh((l  +  i)no)  _  sinh((l  +  0qoy)  x. 

(l  +  /)T|o  (l  +  OTioCOSh((l  +  j)Tlo) 


(3.39) 


D.  THE  PRESSURE  DISTRIBUTION 


Applying  the  boundary  condition  that  v,  =  0  at  >’=0  to  the  continuity  equation 
(Equation  3.32),  an  ordinary  second  order  differential  equation  for  p,(x)  results  as 
follows: 


d^  p,  2— 

— ^  +  m  p.(x)  =  Q  where 
dx 

^2  _ _ (1 + ^)noV^ 

_  ^tanh((l  +  t)no)^j 

(i+j)no 


(3.40) 
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The  boundary  conditions  necessary  for  determining  p,  (x)  are  found  in  the  geometry  of 
the  problem.  As  shown  in  Figure  3. 1,  the  plates  are  flush  with  the  source  of  the  acoustic 
wave  at  x=0.  This  means  that  the  pressure  amplitude  is  at  a  peak  (magnitude  of  one),  and 
the  velocity  amplitude,  and  therefore  the  pressure  gradient,  is  at  a  node  (magnitude  of 
zero).  Thus,  the  appropriate  boundary  conditions  needed  for  determining  the  pressure 

distribution  are  at  x=0,  p,  =iand  Solving  for  p,  (x)  and  applying  these 

boundary  conditions,  the  result  is 

p,(x)  =  icos(mx)  (3.41) 

This  expression  is  then  substituted  back  into  Equations  3.36,  3.37,  3.38,  and  3.39  to  give 
explicit  forms  of  m,  ,  (t)o, Pi ,  and  v,  respectively. 

The  two  cases  studied  in  this  research  are  the  wide  gap  and  arbitrary  gap  cases 
which  differ  bases  on  the  width  of  separation  between  the  plates  as  compared  to  the 
Stokes’s  boundary  layer  thickness  as  discussed  in  Chapter  I.  The  basic  assumption  that  is 
made  in  the  wide  gap  case  is  that,  in  the  pressure  distribution  of  Equation  3.41,  the 
parameter  “m”  is  assumed  to  be  unity  since  m~].  This  simplifies  the  problem 
significantly,  since  the  pressure  distribution  just  varies  as  cos(x)  and  is  independent  of  any 
other  parameter.  For  the  arbitrary  gap  case,  “m”,  is  left  as  a  function  of  T|o,  y,  and  Pr  as 
shown  in  Equation  3.40,  and  the  pressure  distribution  becomes  dependent  upon  those 
three  parameters  as  well  as  on  the  variable,  x. 
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E.  fflGHER  ORDER  BEHAVIOR 


The  higher  order  velocity,  density,  and  pressure  terms,  Mj,  pj,  Vj  andp2.  and 
higher  will  not  be  considered  in  this  study;  however,  the  higher  order  temperature  term, 
<}), ,  is  significant  and  will  be  considered  in  this  section.  The  equation  to  be  solved  in  this 
consideration  is  the  0(1)  energy  equation: 


Energy:  0(1): 


+  u  ^  +  v  ^  =  _L_^ 
dt  '  dx  '  dy  2r\l  Pr  dy^ 


+  2pr„Ec-^«t),^ 


(3.42) 


j 

Again  applying  the  equality,  =  M ^  (— - ),  and  subsequently  neglecting  terms  of 

m 

order  Af^,  Equation  3.42  further  reduces  to 


Energy:  0(1): 


^  +  0  ^  +  M  ^  +  v  ^  =  _i_^  + 
dt  ^  dt  *  9jc  ^  dy  2r\l  Pr  dy^ 


at  dx  T|o  dy 


(3.43) 


Equation  3.43  is  not  complete;  it  is  missing  one  term.  The  extra  term  comes  from 
expanding  the  coefficient  of  thermal  expansion,  P,  into  its  perturbation  form.  By 
definition. 


_ Lr^PJl'k  = _ Lrr^PJ!.'!  ,  Ar*  + 

^  p*  8r* p*  dr*  ”  97*2  “ 


(3.44) 


where: 


AT*  =  T*-T„  =(t)Ar„ 


Utilizing  the  appropriate  scaling  factors  for  p*  and  T*  from  Equation  3.5,  and  the 
corresponding  perturbation  solution  forms  from  Equation  3.1 1,  the  expanded  equation  for 


pr*  becomes 


rmPm 


(3.45) 


Since 


/  r  \  _ 


p.P. 


(  2 )  =  2  for  an  ideal  gas,  and 


_Y-1 

2c  T  ~  ~Y' 

p  nt 


,  Equation  3.45  simplifies 


to 


pj’ =  P.j;{l  -  M(p,  +  +  O(M') +  ...) 


(3.46) 


Substituting  the  expanded  pj*  expression  into  the  energy  perturbation  equation,  an  extra 

0(1)  term  comes  out  of  the  2Ecpr*(— )^  term  from  the  0(M)  portion  of  pr 

M  at 

multiplying  with  (1/Af).  Therefore,  the  extra  0(1)  energy  term  is 
~2P„r„[p,  +  C  .  ^)<|>o  >  making  the  complete  0(1)  energy  equation  as  follows: 

at  2  at 


Energy:  0(1): 

at  at  dx  dy  2^0  Pr  dy^  dt 


(3.47) 


ax  ay 


at 


at 


Substituting  the  relations  from  Equation  3.17  into  Equation  3.47,  yields  the 
following  complex  partial  differential  equation: 
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Energy:  0(1): 


.2il  ,  rz  A  ^2(V  ,  —  ^^0  .2i>  .  —  ^00  _2(f 


Re[2i^^^e  +lp^<^Qe^‘'  +u 


1  3  +v, 

ox 


dy 


-.«1_  >  ^  <1>I0 

2Ti^Pr  dy^ 


+  Re[ 


J  ^  ^12  ^2il  J  ^  p^n  -r  .•=  -2<»  .  Ec  du,  ^2  _2it 


l^lPr  dy^ 
- 


Tlo  oy 


(3.48) 


+  2£cP„7;i7.  -  2p„  -  p„r„/(Y  -  l)(|>„p,.^"  ] 


Two  equations  may  be  extracted  from  Equation  3.48.  One  yields  information  about  the 
first  harmonic  of  (j), ,  which  is  the  real  variable,  (j),^ ;  solution  of  the  other  equation  yields 

the  second  harmonic,  <t),2 ,  which  is  complex.  The  determination  of  0,2  is  not  of  interest 
in  this  study  since  it  does  not  contain  a  steady  temperature  component;  however,  0,q  is 

the  time-averaged  temperature.  Since  it  is  of  great  importance  to  this  research,  the  time- 
averaged  temperature,  0,o ,  is  considered  in  detail  in  the  next  section. 


F.  THE  TIME-AVERAGED  TEMPERATURE,  0io 


In  Equation  3.48,  there  is  one  term  which  obviously  has  no  time  dependence. 


( 


1  a'0 


10 


2q^Pr  ay= 


);  the  remaining  terms  are  multiplied  by  e^".  This  supposed  time 


dependence  would  imply  a  simple  solution  for  0,o  since  it  would  mean  that 


=  0. 


However,  upon  close  inspection  of  the  terms,  it  is  noted  that  most  of  them  have  a 
steady  (time  independent)  term  embedded  in  them.  Taking  the  time  average  of  Equation 
3.48,  the  only  terms  which  drop  out  (average  to  zero)  are  2i0,2e^", 
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however,  the  calculus  and  algebraic  manipulations  involved  in  solving  the  equation  by 
hand  are  extremely  cumbersome.  For  example,  evaluating  the  time  average  of  the  first 
term,  <  ip  ,<!)(,  > ,  is  accomplished  using  the  definition  of  time  average  from  Equation  3.50 


and  substituting  into  it  the  expressions  for  p,  and  (l)^  from  Equations  3.37,  3.38,  and 
3.41  as  follows: 
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>=  \ Re[p, H)(|).  »  ♦]  =  1  Re(icos(/>a)[l  +  ♦ 

2  2  cosh((l  +  j)noVPr) 


(-0[2/cos(mx)(l  - 


cosh((l  +  i)^\oy^/PT) 
cosh((I  + /)r|o  VPt) 


)]**} 


where  m  = 


_ 0  +  i^ovPr 


^tanh((l  +  /)no)^j 

0  +  i)no 


(3.51) 


Once  the  complex  conjugate  of  <|)o  is  correctly  computed,  the  real  part  of  the  entire 
expression  must  be  evaluated  without  error.  This  establishes  the  first  time  average  term, 
but  there  are  six  more  which  must  be  obtained  in  a  similar  fashion.  The  summation  of  all 
seven  resulting  time  average  terms  must  then  be  integrated  twice,  applying  appropriate 
boundary  conditions,  to  obtain  an  expression  for  The  vast  amount  of  careful 

mathematical  manipulation  and  computation  required  to  solve  for  ((),o  points  directly  to 

the  use  of  a  powerful  symbolic  mathematical  software  tool.  Thus,  MAPLE,  just  such  a 
computer  software  package,  is  utilized  to  solve  for  (|>,o ,  given  all  of  the  other  expressions 

for  M, ,  <t)o ,  p , ,  V, ,  and  p, .  (j),^  is  found  first  for  the  simpler  wide  gap  case  and  then  for 
the  arbitrary  gap  case.  The  resulting  expressions  for  (t),o  are  much  too  lengthy  to  include 
explicitly  in  this  paper;  however,  the  temperature  distributions  described  by  <j>,o  are 
plotted  and  discussed  in  Chapter  IV. 
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IV.  RESULTS  AND  DISCUSSION 


The  results  of  this  work  are  comprised  of  the  solutions  to  the  governing  equations 
as  given  in  Chapter  HI  along  with  plots  of  these  solutions.  Tliese  plots  are  presented  and 


discussed  in  the  following 


thermal  penetration  depth  is 


sections  for  different  values  of  Tio>/ft=—  where  the 

5„ 


2oc  S 

5„  = .  —  =  .  The  parameter,  r\.^fPr  ,  has  been  found 

V  CO  VPr 


to  be  the  more  appropriate  gap  width  parameter  than  T|o  in  dealing  with  thermal  effects  in 
such  oscillatory  flow  problems.  Since  the  form  of  the  leading  order  pressure  distribution, 
Pi(x),  is  crucial  in  determining  the  behavior  of  the  remaining  leading  order  solutions  and 
the  higher  order  solutions,  the  pressure  distribution  results  are  analyzed  first.  Then,  the 
leading  order  velocity,  temperature,  and  density  solutions  and  the  higher  carder  time- 
averaged  temperature  term  are  presented.  Finally,  the  recovery  factor,  suitably  defined 
for  this  model  problem,  is  presented,  plotted  and  discussed. 


A.  THE  PRESSURE  DISTRIBUTION 


The  pressure  distribution,  Pi(x),  for  this  model  problem  is  found,  as  given  in 
Chapter  ID,  as 


Pi(x,t)  =  Re[p,(A:)e"],  where  /?,(x)  =  icos(mx)  and 


m  = 


_ (1  ^)T|q  yPr 


^tanh((l-n-)no)^j 

0+»)no 


(4.1) 
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Note  that  the  eigenvalue  “m”  is  in  general  complex;  however,  with  increasing  gap  width, 
the  parameter  “m”  approaches  a  value  of  one.  Therefore,  in  the  idealized  wide  gap  case, 
it  is  assumed  that  m=l,  and  the  pressure  distribution  is  a  simple  cosine  function,  i.e. 
p,(x)  =  icosx  and  Py{x,t)  =  -cosxsint .  In  general,  for  the  arbitrary  gap  case,  “m”  is  a 

complex  number,  m  =  a  -  bi.  As  gap  width  (or  iio)  increases,  “a”  reduces  to  a  value  of 
one  and  “fe”  to  a  value  of  zero,  approaching  the  results  of  the  idealized  wide  gap  case. 
For  all  of  the  trials  calculated,  the  values  of  “m”  are  given  in  Table  4. 1 .  The  value  of  “m” 
also  depends  on  the  gas  properties,  y  and  Pr.  Thus,  for  the  same  rjo  value  or  gap  width, 
the  magnitude  of  “m”  will  vary  for  different  gases.  The  ideal  gas  utilized  in  this  work  is 
air;  therefore,  the  property  values  used  are  7^1.4  and  Pr=0.7. 


Table  A 

kl  Values  of  “m”  in  Pressure  Distribution. 

1 

o 

m 

M 

1 

1.3047-0.5939/ 

1.443 

1.5 

1.1814-0.3165/ 

1.223 

2 

1.1405-0.2055/ 

1.159 

3 

1.1002-0.1229/ 

1.107 

4 

1.0762-0.0884/ 

1.080 

5 

1.0613-0.0688/ 

1.064 

Two  plots  of  the  pressure  distribution,  pi(x,f),  are  given  in  Figure  4.1,  for  iloVft 
values  of  1.5  and  5.  Considering  the  Tj^VPir  =5  case,  the  pressure  distribution  closely 
resembles  the  idealized  wide  gap  case  of  p^{x,t)  =  -cosxsinr .  As  gap  width  decreases, 
the  amplitude  of  the  cosine  function  increases,  and  the  waveform  is  considerably  distorted 
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as  illustrated  by  the  tIqVPt  =1.5  plot.  This  variation  in  amplitude  is  clearly  observed  in 

Figure  4.2  where  Pi{x,t)  is  plotted  at  two  particular  times  t=—  and  t=n  for  ilj,  Vft  values 

of  1.5,  2,  3  and  5.  To  emphasize  the  distortion  as  compared  to  the  idealized  wide  gap 
case,  the  function,  p^{x,t)  =  -cosxsint  is  also  plotted  in  Figure  4.2.  Based  on  the 
calculated  values  of  “m”  and  the  plots  of  the  pressure  distributions,  the  gap  width 
corresponding  to  tjoVPr  =5  is  determined  to  the  sufficiently  close  to  the  idealized  wide 
gap  case  to  be  considered  as  a  “wide  gap”.  Again  as  r|o  decreases  into  the  “narrow  gap” 
region  (rioVPr  <5),  the  pressure  distribution  significantly  deviates  from  the  simple  cos{x) 
behavior  of  the  idealized  wide  gap  case. 

B.  LEADING  ORDER  BEHAVIOR 


The  leading  order  first  harmonic  solutions  for  the  primary  variables  (velocity 
components,  temperature  and  density)  are  given  as  follows  (as  discussed  in  Chapter  HI): 

M,  (x,  y,  t)  =  Re[M;  {x,  y)e'‘  ] ,  <|)o (j:,  y,  0  =  Re[0o ix,  y)c" ] ,  p ,  (jc,  y,  r)  =  Re[p ,  (x,  y)e"  ] ,  and 


Vi(A:,y,0  =  Re[v,(x,y)e"],  where 


dx  cosh((l  +  j)no) 


(t)o(j:,y)  =  2p,(x)[l- 


cosh((l  +  i)noyV^) 

cosh((l  +  /)TioVPr)  ^ 


PiU.y)  =  p,(j^)[Y-(Y-i)(i- 


cosh((l  +  i)noyVPr) 
cosh((l  +  i)noVPr) 


(4.2) 

(4.3) 


(4.4) 
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V,  {X,  y)  =  ip,  (jf)[(  1  -  3')  +  (y  - 1)( 


tanh((l  +  /)noVi^) 

(I  +  OiloV^ 


- s±nm^iyn,y^)  (45) 

(l  +  i)T|(,vftcosh((l  +  i)T|oVFt)  dx^ 

^  ^tanh((l+OTio)  _  sinh((l  +  OTioy) 

(l  +  /)ilo  (l  +  l)T|oCOSh((l  +  OTlo)'^'^ 

These  quantities,  as  opposed  to  the  pressure  distribution,  depend  on  both  spatial 
variables,  x  and  y.  Their  x  dependence  is  given  by  either  cos(mx)  or  sin(mx)  which  steins 
from  their  dependence  on  the  pressure  distribution  or  the  pressure  gradient.  Because  of 
this  periodicity  in  x,  all  of  the  information  about  the  x  dependence  of  a  particular  leading 
order  distribution  is  contained  in  an  interval  of  one  half  of  the  wavelength  (i.e. 
0<x<7t ).  In  the  j  direction,  the  quantities  <|>o,  and  p,  are  symmetric  about  the 
centerline  of  the  channel  (y=0),  whereas  the  velocity  in  the  y  direction,  v,,  is 
antisymmetric  about  the  centerline.  Therefore,  the  results  for  all  of  these  distributions  are 
determined  and  presented  for  one  half  of  the  channel  only,  specifically  for  0  <  y  <  1 .  The 

plots  of  the  leading  order  density  distribution,  pi(x,y,0,  are  not  included  in  the  discussion 
since  the  behavior  of  the  density  function  is  not  of  direct  interest  to  this  work. 

The  velocity  distribution,  u,(x,y,t),  is  given  in  Figure  4.3  for  ;c=JC  and  TioVPr 

values  of  1.5  and  5.  Ui(x,y,t)  depends  on  the  pressure  gradient;  therefore,  it  varies  in  a:  as 
sin(/nx).  This  distribution,  as  observed  in  Equation  4.2,  varies  in  y  as 

Re[l - .  Figure  4.4  shows  Ui(x,y,t)  for  x=ii  at  two  instants  in  time,  t=— 

cosh((l  +  i)no)  2 

and  t=K  for  TioVPr  =  1.5,  2,  3,  and  5.  The  distributions  clearly  show  that  the  velocity 


vanishes  at  the  wall  (y=l)  and  thereby  satisfies  the  no  slip  condition.  Also,  because 
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Ui{x,y,t)  is  symmetric  about  the  centerline  of  the  channel,  the  gradient  there  is  equal  to 
zero  as  shown  by  the  zero  slope  at  y=0.  For  wide  gaps  or  large  values  of  t|o,  the  velocity 
component,  u,,  increases  in  magnitude  from  zero  at  the  wall  (y=l)  to  a  maximum  value 
and  then  decreases  slightly  to  level  off  with  a  slope  of  zero  at  the  centerline  Cy=0).  As  the 
gap  narrows  or  tio  decreases,  the  m,  velocity  value  increases,  and  the  location  of  the  peak 
velocity  in  the  channel  shifts  toward  the  centerline.  In  fact,  for  rioVft=1.5,  the 

maximum  velocity  value  occurs  at  the  centerline,  and  the  distribution  is  more  parabolic  in 
shape  (similar  to  that  of  Poiseuille  flow).  These  trends  are  clearly  observed  in  the  curves 

7C 

for  both  snapshots  in  time,  t=—  and  r=7t,  in  Figure  4.4. 

The  leading  order  temperature  solution,  (|)o(jc,y,0,  is  illustrated  in  Figure  4.5  for 
x=n  and  =  15  and  5.  The  x  dependence  is  a  cosine  function,  jcos(mx),  which  is 

caused  by  the  direct  relationship  of  ({)o  to  the  pressure  distribution.  The  y  dependence  of 
^  cosh((l +  /)Ti„yVft)^, 

<|)o  IS  of  the  form  Re[i(l - — — .  ”  /^  )]  as  observed  in  Equation  4.3.  In  Figure 

cosh((l  +  i)r|oVPr)  ^  ® 

4.6,  where  (|)o  is  plotted  for  x=Tt,  at  two  instants  in  time,  t=—  and  t=K,  and  for  tIqVPt 

values  of  1.5,  2,  3,  and  5,  the  dependence  of  (jio  on  y  and  rjo  is  displayed,  and  its  similarity 
to  that  of  the  Ui  velocity  distribution  is  noted.  The  boundary  conditions  of  no  slip  at  the 
wall  and  symmetry  at  the  centerline  are  clearly  met  as  <j)o=0  at  y=l  and  its  slope  vanishes 
at  y=0.  Considering  the  curves  for  time  t=it,  for  large  gap  widths  or  Tjo  values,  the 
temperature,  <t>o,  increases  in  magnitude  from  zero  at  the  wall  to  a  peak  and  then  decreases 
to  the  centerline  temperature  value  which  is  approached  at  zero  slope.  Again,  as  the  gap 
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width  or  Tjo  value  decreases,  the  peak  shifts  toward  the  centerline  and  its  magnitude 
increases.  As  with  the  Uj  velocity  component,  the  shape  of  this  temperature  distribution 
becomes  more  parabolic  as  the  gap  narrows  or  T|o  decreases.  The  dependence  of  the 

K 

temperature,  (t)o,  on  gap  width  or  Tjo  is  not  accurately  portrayed  in  the  plots  for  time, 

This  is  explained  by  the  nature  of  the  time  dependence  of  the  temperature  distribution, 
^o(x,y,t).  A  phase  lag  occurs  in  time,  and  the  amount  of  this  phase  lag  depends  on  the 

value  of  rio;  therefore,  the  temperature  distributions  at  the  instant  of  time,  t=—,  for  the 

different  gap  widths  do  not  correspond  to  the  same  time  positions  in  their  respective 
periods. 

The  velocity  component  in  the  y  direction,  \i(x,y,t),  is  presented  in  Figure  4.7  for 
x=n  and  T|oV^  =  1.5  and  5.  From  Equation  4.5,  it  is  observed  that  Vi(x,y.r)  depends  on 

both  the  pressure  distribution  and  its  second  derivative.  Therefore,  the  x  dependence  is  a 
cosine  function,  cos(mx).  The  y  dependence  of  V|(x,y,0  is  a  combination  of  linear  terms, 
hyperbolic  sine  terms  and  constants,  resulting  in  a  function  which  is  antisymmetric  about 
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y=0.  In  Figure  4.8,  v,(x,y,0  is  plotted  for  jt=7t,  at  two  instants  in  time,  t=—  and  t=n,  and 

2 


for  tIqVPt  values  of  1.5,  2,  3,  and  5.  The  boundary  conditions  of  no  penetration  at  the 
wall  and  “zero  crossover”  at  the  centerline  for  an  antisymmetric  distribution  are  satisfied 
as  v,=0  at  y=l  and  y=0  respectively.  The  behavior  change  in  V|  as  gap  width  or  t|o 
changes  is  clearly  illustrated  by  the  t=Tt  curves.  For  large  gap  widths,  the  y  component  of 
the  velocity,  v,,  increases  in  magnitude  from  zero  at  the  centerline  to  a  peak  near  the  wall 
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and  then  decreases  rapidly  in  magnitude  back  to  zero  at  the  wall.  An  inflection  point  is 
observed  between  the  centerline  and  the  peak  value  for  wide  gaps.  As  the  gap  narrows  or 
Tio  decreases,  the  peak  increases  in  magnitude  and  shifts  away  from  the  wall.  Also,  the 
shape  of  the  distribution  becomes  more  parabolic  for  narrower  gaps  (i.e.  the  inflection 
point  observed  in  the  wide  gap  trials  disappears  as  Tio  decreases).  The  trend  cannot  be 

K 

distinguished  in  the  t=—  plots  for  reasons  relating  to  the  nature  of  the  time  dependence  of 
v,(jc,y,r)  which  are  similar  to  those  discussed  in  the  (t>o(^,y,0  results. 

C.  THE  TIME-AVERAGED  TEMPERATURE 


The  time-averaged  temperature  distribution,  (|)io(Jc,y),  is  shown  in  Figure  4.9  for 
=  1.5  and  5.  From  the  partial  differential  equation  for  <|)io  as  presented  in  Chapter 
in  and  again  here  for  convenience. 


<  ip,<l)o  >-!-<«,  >  +  <  V,  > 


I  9'(t) 


10  ,  ^  /-^“l  \2 


9jc 

TT  ^1 


+  <• 


ay  2Ti^Pr  r\l  dy 


■(^r  > 


(4.6) 


+  <  -^>-<  2P„7;/p,p,  >  -  <  p„r„i(y  -  i)(t)oP,  > 


the  X  dependence  is  determined  by  the  various  time-averaged  terms  in  the  forcing 
function.  The  terms  involve  the  multiplication  of  two  of  the  leading  order  first  harmonic 
solutions  of  the  primary  variables,  <()o,  p,  ,  and  v,,  or  their  derivatives.  Substituting  the 
appropriate  distributions  from  Equations  4.2,  4.3,  4.4,  and  4.5  into  Equation  4.6,  the  x 
dependence  of  <l))o  is  observed  to  stem  from  terms  where  the  pressure  distribution  is 
multiplied  by  itself  as  well  as  terms  where  the  pressure  gradient  is  multiplied  by  itself. 
Therefore,  the  resulting  time-averaged  temperature  distribution  varies  in  jr  as  a 
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combination  of  sin^(mx)  and  cos^(mx)  terms.  In  Figure  4.10,  the  time-averaged 

temperature,  (1),o(jc,j),  is  shown  for  x=0,  x=—  and  x=7i  and  for  'noVft  values  of  1.5,  2,  3, 

2 

and  5.  The  temperature,  ^lo,  vanishes  at  the  wall  Cy=l),  satisfying  the  no  slip  boundary 
condition;  and  its  gradient, =  0  at  the  centerline  (y=0),  thereby  satisfying  the 

symmetry  condition.  From  the  jc=0  curves,  the  variation  in  <(>io  based  on  T|o  is  clearly 
observed.  For  wide  gaps  or  large  Tio  values,  the  time-averaged  temperature  increases  from 
zero  at  the  wall  to  a  maximum  value  fairly  close  to  the  wall  and  then  decreases  gradually, 
leveling  off  to  a  zero  slope  at  the  centerline.  Thus,  for  wide  gaps  most  of  the  temperature 
variation  occurs  near  the  wall.  As  the  gap  size  becomes  smaller  or  tIo  decreases,  the 
magnitude  of  the  time-averaged  temperature  gradient  decreases,  and  the  temperature 
variation  is  spread  more  and  more  evenly  across  the  channel  width.  Further  discussion  of 
the  affect  on  the  time-averaged  temperature  gradient  of  varying  the  gap  width  is  presented 
in  the  next  section  on  the  recovery  factor.  The  time-averaged  temperature  distribution, 
<|>io(Jc,y),  is  not  symmetric  in  x.  This  is  illustrated  by  the  fact  that  the  plots  for  jc=7C  do  not 
overlap  with  the  curves  for  Ar=0. 

D.  THE  RECOVERY  FACTOR 

More  important  than  the  results  for  the  time-averaged  temperature  distribution, 
<t>io,  presented  above,  it  is  the  gradient  of  this  temperature  at  the  wall  which  allows  key 
deductions  to  be  made  about  the  time-averaged  heat  transfer  effects  in  this  problem.  In 
order  to  interpret  these  heat  transfer  results  in  a  more  conventional  manner,  they  are 
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presented  in  terms  of  a  suitably  defined  recovery  factor.  The  recovery  factor,  r,  is  the 
ratio  of  recovered  thermal  energy  at  an  adiabatic  surface  to  the  fluid’s  kinetic  energy  as 
discussed  in  Chapter  II.  The  form  of  “r”  utilized  in  this  research  conforms  to  convention 
and  is  given  by 

^0 


T„  is  the  fixed  temperature  of  the  upper  plate.  Taw  is  the  adiabatic  wall  temperature  and  it 
is  equal  to  the  value  of  T*  at  y*=-yo  when  the  wall  there  is  adiabatic.  The  adiabatic 

dT* 

condition  is  that  the  temperature  gradient,  - —  =  0  at  the  wall  in  question,  y*=-yo.  In  the 

ay* 

present  problem.  Taw  would  correspond  to  the  time-averaged  temperature  achieved  by  a 
wall  if  it  were  made  adiabatic.  Let  the  corresponding  dimensionless  temperature  for  this 
adiabatic  wall  be  denoted  in  general  as  tj/.  Then  it  is  the  corresponding  time-averaged 
temperature,  ,  at  the  adiabatic  wall  that  is  of  interest  here.  Expressing  the  recovery 

factor  and  the  adiabatic  conditions  in  terms  of  the  dimensionless  parameters  yields: 


at  y  =  -l, 

r^=\j/,oAr„^-hr„  at  y  =  -l, 
c^(x|/,oAr„,+7;-rj  i|t.„Ar„ 


r-^=Vio» 


or  r=\|/,o(y  =  -l) 

where  \jr  is  the  dimensionless  time-averaged  temperature  obtained  in  the  adiabatic  wall 


47 


In  order  to  find  “r”  in  this  manner,  however,  it  is  necessary  to  additionally  solve 
for  the  temperature  distribution,  m ,  in  the  adiabatic  wall  case.  To  avoid  this  extra  time- 

consuming  step,  a  relationship  will  be  proven  between  the  value  of  the  time-averaged 
temperature,  ,o ,  at  an  adiabatic  wall  and  the  value  of  the  time-averaged  temperature 


gradient,  the  wall  in  the  fixed  temperature  case  that  is  the  subject  of  this 


research.  The  two  distributions  appear  to  have  approximate  shapes  as  given  in  Figure 
4.11. 


Figure  4. 1 1  Approximate  Temperature  Distributions  for  Two  Cases: 
(a)  Fixed  Wall  Temperature  and  (b)  Adiabatic  Wall  Temperature 


T  * 

<t>.  - — ^ 


■“  Yi 

Boundary  Conditions: 

<|),o=0  at  y  =  ±l 


Vio  = 


y  ♦  — y 


2c. 


'10 


dy 


0  at  y  =  0 


Boundary  Conditions: 
\|r,o=0  at  y  =  l 
5¥io 


dy 


=  0  at  y  =  -1 
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The  recovery  factor  has  already  been  shown  in  Equation  4.8  to  be  equal  to 

V io(>’  =  -1)  •  It  will  now  be  proven  that  V|/  ,o(y  =  -1)  =  2^^(y  =  -1) .  From  the  nature 

dy 

of  the  partial  differential  equation  for  0,0  in  Equation  4.6,  it  is  known  that  the  y 
dependence  of  0,o  can  be  expressed  as 

<l>io(>')  =  /(>’)  +  Ci>'  +  C2  (4.9) 

where  from  the  nature  of  the  forcing  function  in  Equation  4.6,  it  can  be  shown  thaty(y)  is 
an  even  function.  Now,  applying  the  appropriate  boundary  conditions  at  the  walls  for  the 
present  problem  yields 

‘l>.oa)  =  0  =  /a)  +  C,+C2  and 
0,o(-l)  =  O  =  /(-l)-C,+C2 

Since is  an  even  function, /(I )=/(-!),  and  the  constants  of  integration  are  determined 

to  be  C,  =  0  and  Cj  =  -/(I) .  Thus,  the  resulting  expressions  for  0,o  and  are 

dy 

<l>io(>')  =  /(y)-/(l)  (4.11) 

^  =  /'W  (4.12) 

Similarly,  for  the  adiabatic  wall  problem,  the  time-averaged  temperature  has  the 
same  governing  equation  since  the  primary  leading  order  variables  are  the  same.  Then, 
\|/  ,0  is  expressed  as  the  same  function  of  y  but  with  different  constants  of  integration  as 
follows: 

^lo(y)  =  fiy)  +  D^y  +  D2  (4.13) 
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Applying  the  appropriate  boundary  conditions  for  this  case,  which  are  y  =0  at 


y=l  and  =0  aty=-l,  yields 
By 


V,o(l)  =  0  =  /(l)  +  D,  +  D2  and 
^(1)  =  0  =  /'(1)  +  D, 


(4.14) 


Thus,  the  constants  are  found  to  be  Z),  = -/'(-I)  and  Dj  = /'(-1)-/(1),  and  the 
resulting  function  is  found  to  be 


¥  ,o(>’)  =  fiy)  -  /'(-1)3' + /'(-I)  -  /(I) 


(4.15) 


Evaluating  the  adiabatic  wall  temperature,  \|f  ,0  aty=-l,  yields 


¥  ,o(-l)  =  /(-I)  +  /'(-I)  +  /'(-I)  -  /(I) 


(4.16) 


Again,  noting  thaty(y)  is  an  even  function  and  using  the  properties  of  even  functions,  the 
previous  equation  is  reduced  to 


\|r,o(-l)  =  2/'(-l) 


(4.17) 


Also,  from  Equation  4. 12,  it  is  seen  that  -^(y  =  -1)  =  /'(-I) .  Therefore, 

By 


xir,o(-l)  =  r  =  2^(y  =  -l) 
By 


(4.18) 


This  proof  allows  “r”  to  be  defined  and  calculated  in  terms  of  the  temperature 
distribution,  <|),o ,  which  has  been  determined  in  the  present  problem  without  having  to 


redo  all  of  the  work  necessary  to  solve  the  adiabatic  wall  case.  The  gradient,  ,  is 

By 

easily  obtained  and  evaluated  at  the  wall  and  then  doubled.  Since /(y)  is  an  even  function, 
/'(y)is  an  odd  function;  therefore,  it  is  known  that/ '(-1)  =  -/'(!).  The  domain  of 
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analysis  is  0  <  y  <  1 ;  therefore,  for  convenience,  the  recovery  factor  calculated  in  this 
research  is  determined  by 

/•  =  -2^(y  =  l)  (4.19) 

oy 

The  resulting  values  of  “r”  may  be  positive  or  negative.  Positive  values  of  ‘V’ 
correspond  with  the  physical  situation  of  the  wall  being  heated  by  the  fluid.  ITiis 
phenomenon  of  the  walls  being  heated  by  viscous  dissipation  is  observed  in  high  speed 
mean  flows  such  as  in  the  example  cases  presented  in  Chapter  II.  Negative  values  of  the 
recovery  factor  indicate  that  the  direction  of  heat  flow  is  away  from  the  wall  as  shown  in 
Equation  4.19.  In  other  words,  the  wall  is  actually  being  cooled  by  the  fluid  in  dtose 
instances.  This  cooling  potential  is  a  very  interesting  attribute  of  high  speed  oscillatory 
flows,  and  it  is  this  phenomenon  which  is  the  underlying  principle  that  makes 
thermoacoustic  refrigeration  possible.  Thermoacoustic  heating  and  cooling  effects  caused 
by  thermoacoustic  streaming  in  a  resonance  tube  were  predicted  theoretically  and  tfien 
validated  experimentally  by  Merkli  and  Thomann  (1975).  A  subsequent  review  of  die 
topic  was  provided  by  Rott  (1980).  Experimental  results  were  also  obtained  by  Wheatley, 
et  al.  (1983)  which  verified  the  existence  of  this  cooling  effect  in  a  stack  of  plates  placed 
in  an  acoustically  resonant  tube. 

The  recovery  factor,  r,  is  calculated  using  Equation  4.19  and  then  plotted  for  each 
trial  run  for  this  problem.  The  results  are  illustrated  in  Figure  4.12  for  one  half  of  a 
wavelength,  0  <  x  <  7t ,  and  for  Tip  VPr  =1 .5,  2,  3,  and  5.  For  each  gap  width  examined, 
the  recovery  factor  takes  on  both  positive  and  negative  values  which  corresponds  to 
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heating  and  cooling  of  the  channel  wall  respectively.  In  general,  at  jc=0,  the  recovery 
factor  is  a  positive  local  maximum.  As  jc  increases,  “r”  decreases  for  each  gap  width  until 
it  reaches  a  minimum  value  which  is  negative.  The  recovery  factor  then  increases  as  x 
increases  until  it  reaches  its  maximum  positive  value  before  again  beginning  to  decrease. 

As  observed  in  Figure  4.12,  the  recovery  factor  is  generally  higher  in  magnitude 
for  all  X  values  plotted  and  the  curve  is  steeper  (i.e.  the  gradients  are  higher  in  magnitude) 
for  a  wide  gap  than  for  a  narrow  gap.  The  fact  that  the  “r”  values  are  higher  for  wide 
gaps  than  for  narrow  gaps  corresponds  with  the  higher  time-averaged  temperature 
gradients  observed  near  the  wall  (y=l)  for  wide  gaps  in  Figure  4.10.  According  to 
Fourier’s  law  as  well  as  the  meaning  of  the  recovery  factor  as  described  earlier,  these 
observations  would  imply  greater  heat  fluxes  for  greater  gap  widths.  This  apparent 
implication  does  not  agree,  however,  with  what  has  been  experimentally  observed. 
Therefore,  in  order  to  better  understand  physically  what  “r”  is  measuring  in  this  model 
problem  of  oscillatory  flow,  the  equation  for  the  recovery  factor  utilized  in  this  work. 
Equation  4.19,  is  expanded  by  substituting  back  into  it  its  dimensional  quantities.  From 
Equation  4.19,  the  recovery  factor  is  proportional  to  the  dimensionless  time-averaged 
temperature  gradient  as  follows: 


r  oc 


By 


(4.20) 


Using  the  appropriate  scaling  relations  involving  (|)|„  and  y  from  Equations  3.4  and  3.5, 
the  relationship  in  Equation  4.20  becomes 
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(4.21) 


yp  ^^0  * 

X.  ^y* 

K 

Manipulating  the  expression  in  Equation  4.21  in  order  to  relate  the  recovery  factor  to  the 
physical  parameters  which  are  more  appropriate  for  this  oscillatory  flow  problem 
produces 


4)(f ) 

c  2c„ 


(4.22) 


where  ■^  =  'noVPr  is  again  the  gap  width  parameter  utilized  in  this  work.  Substituting 
Oa 

y 

the  relations,  M  =  —  and  =yRT^  ,  into  Equation  4.22,  the  expression  for  the 
c 

recovery  factor  becomes 


a(%^  TioVpr  dC%) 

M4^)  di%) 

2  2 


(4.23) 


From  Equation  4.23,  the  apparent  contradiction  to  physics  observed  in  Figure  4.12 
is  explained.  It  is  clear  from  Equation  4.23  that  for  increasing  values  of  TioV^ ,  an 

increase  in  the  magnitude  of  “r”  does  not  automatically  imply  an  increase  in  the  time- 
averaged  temperature  gradient  on  the  right  hand  side  (RHS).  In  order  to  distinguish 
clearly  the  effect  on  the  time-averaged  temperature  gradient  as  the  gap  width  (or  TipVPr ) 
changes,  the  recovery  factor,  r,  must  be  scaled  by  T|oVft  as  follows: 
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r 


(4.24) 


_  1  ^C'%) 

n.VpF  “  3(’X) 

2 

Therefore,  for  a  given  Mach  number,  M,  and  for  a  particular  gas  (i.e.  y  fixed),  the 
r 

quantity, - ,  in  Equation  4.24,  carries  the  physical  information  about  how  the  time- 

TloVPr 

averaged  temperature  gradient  (and  thus  the  heat  flux)  varies  with  different  gap  widths  (or 
values  of  tipV^ ). 

The  recovery  factor,  scaled  by  TloV^  ,  is  plotted  in  Figure  4.13  for  =  1-5, 

2,  3,  and  5  and  illustrates  the  true  heating  and  cooling  effects  at  the  wall  (y=l).  For  a 

wide  gap  as  compared  to  a  narrow  gap,  the  “scaled  recovery  factor”  ( — ^7=)  is  lower  in 

magnitude  when  x=0  and  remains  lower  until  approximately  x=  0.27c.  At  this  point  all  of 
the  curves  cross  so  that  the  scaled  recovery  factor  is  slightly  higher  in  magnitude  for  wide 

gaps  until  it  crosses  zero.  During  the  initial  cooling  portion  of  the  curve  ( — ^=<0),  the 

TloVPr 

wide  gap  magnitudes  are  lower  again  until  the  curves  turn  at  their  respective  minima. 
Once  the  heating  region  is  reached  again  ( — ^7=>0)  for  all  of  the  curves,  the  scaled 

recovery  factor  is  again  lower  in  magnitude  for  wide  gaps.  The  — curves  are  flatter 

TioV^ 

(i.e.  gradients  are  lower  in  magnitude)  for  wide  gaps  for  generally  all  values  of - 7= . 

Considering  the  portion  of  the  curves  where  the  scaled  recovery  factor  is  negative,  the 
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following  observations  are  noted.  For  ri(,VPr=5,  the  zero  crossings  of  the  scaled 

recovery  factor  occur  at  approximately  x=0.36ic  and  jc=0.587C.  Thus,  the  length  of  the  x- 
r 

interval  for  which  -  is  negative  is  approximately  0.22jc.  As  the  gap  width  narrows 

or  Tjo  decreases,  the  average  magnitude  of  the  scaled  recovery  factor  decreases  and  the 
magnitude  of  its  gradient  increases.  This  results  in  a  shift  of  both  zero  crossings  to  lower 
X  values  as  rio  decreases  as  well  as  a  shortening  of  the  x-interval  over  which  the  scaled 
recovery  factor  is  negative.  For  example,  when  T|oVPr  =1.5,  the  zero  crossings  of  the 

scaled  recovery  factor  occur  at  approximately  x=0.33n  and  x=0.467i,  resulting  in  an  x- 
r 

interval  where  the  ^  value  is  negative  of  approximately  0.1371. 

The  higher  magnitudes  of  the  positive  scaled  recovery  factor  for  narrower  gaps 
translate  into  higher  heat  fluxes  at  the  wall  of  the  channel  for  those  x  locations.  Just  as  in 
the  previous  discussion  concerning  the  conventional  recovery  factor,  r,  for  the  x-intervals 
where  the  scaled  recovery  factor  is  positive,  the  channel  wall  is  being  heated  by  the  fluid 

through  viscous  dissipation.  For  the  x-interval  where  — ^7=  is  negative,  the  channel 

wall  is  actually  being  cooled  by  the  oscillating  fluid.  The  gradient  of  — ^-r=  is  also 

tIoVpF 

important  because  it  determines  the  heat  pumping  capacity  available.  More  specifically, 

for  the  same  magnitude  of - y=  (or  cooling  effect),  a  higher  gradient  equates  to  more 

TioVPr 

heat  pumping  capacity  available  for  the  same  length  of  plates. 
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This  cooling  effect,  produced  in  zero-mean  oscillatory  flows,  is  utilized  in 
thermoacoustic  refrigeration  applications.  The  stack  of  the  thermoacoustic  refrigerator  is 
positioned  in  a  tube  where  part  of  the  stack  is  in  a  region  of  x  values  where  the  scaled 
recovery  factor  is  positive,  and  the  remaining  portion  of  the  stack  is  in  a  region  where 

- IS  negative.  The  end  where  cooling  occurs  (negative  - p=  values)  is  where 

the  refrigerating  effect  is  produced.  The  heat  generated  at  the  other  end  of  the  stack 
(where  the  scaled  recovery  factor  is  positive)  is  removed  via  a  heat  sink.  The  stack  region 
of  a  typical  thermoacoustic  refrigerator  is  very  short  compared  to  the  wavelength  of  the 
acoustic  wave  producing  the  oscillatory  flow;  therefore,  an  analysis  of  the  scaled  recovery 


factor  in  the  region  of  x  where 


r 

tIoVpF 


just  becomes  negative  is  important  for  this 


application. 


Comparing  both  the  magnitudes  of - 

TloVPr 


and  the  gradients  in  the  region 


of  X  from  approximately  0.33n  to  approximately  0.47t  in  Figure  4.13,  it  is  observed  that 
the  highest  magnitude  and  gradient  of  the  scaled  recovery  factor  occur  for  the  case  where 
TjoVPr  =2.  Therefore,  for  a  stack  place  with  its  cooling  end  in  that  x  region,  the  heat  flux 


away  from  the  wall  and  the  heat  pumping  capacity  are  both  greater  for  a  gap  width 
corresponding  to  T|q  VPt  =2  than  for  any  of  the  other  gap  widths  examined. 
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V.  CONCLUSION 


This  work  involved  finding  the  analytical  solution  to  a  fundamental  problem  in 
order  to  study  the  time-averaged  heat  transfer  effects  induced  by  high  speed  zero-mean 
internal  oscillatory  flows.  The  governing  equations  of  continuity,  momentum,  energy, 
and  state  were  established  for  the  model  problem  which  consisted  of  an  ideal-gas  host 
fluid  between  two  parallel  plates  excited  into  a  zero-mean  oscillatory  flow  by  a  high- 
intensity  standing  acoustic  wave. 

Introducing  the  proper  scaling  relations,  the  governing  equations  were 
transformed  into  dimensionless  equations;  then,  order  of  magnitude  analyses  were 
conducted  to  simplify  the  resulting  equations.  Employing  perturbation  solution  methods 
with  the  Mach  number,  M,  used  as  a  perturbation  parameter,  the  equations  were  expanded 
and  separated  based  on  the  orders  of  M.  The  most  significant  order  equations  were  then 
analytically  solved  using  complex  variable  methods  for  the  unknown  velocity,  density, 
temperature,  and  pressure  terms.  The  final  aim  was  to  determine  the  time-averaged 
temperature  distribution  for  various  gap  widths  in  order  to  analyze  the  heat  transfer 
effects  induced  by  the  zero-mean  oscillatory  flow. 

The  analysis  was  completed  using  air  as  the  host  fluid  for  gap  widths 
corresponding  to  various  TjoVft-  values  up  to  T|oVPr  =5,  which  was  determined  to  be  the 

lower  limit  for  consideration  as  a  “wide  gap”.  The  leading  order  primary  variables  were 
presented  graphically  and  discussed  in  detail  in  Chapter  IV  along  with  the  time-averaged 
temperature  and  corresponding  recovery  factors.  According  to  the  interpretation  of  these 
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results,  using  air  as  the  host  fluid,  the  most  suitable  choice  of  gap  width  in  terms  of  the 
cooling  effects  applicable  to  thermoacoustic  refrigeration  corresponded  to  a  value  of 
TloVPr=2. 

This  was  a  fundamental  study,  hoped  to  be  a  foundation  upon  which  to  be  built,  in 
the  area  of  recovery  factors  in  zero-mean  internal  oscillatory  flows.  Further  work  is 
suggested  such  as  determining  the  most  suitable  gas  to  be  used  as  a  host  fluid  (i.e.  to 
maximize  the  cooling  effects  produced).  Also,  model  problems  involving  different 
geometries  (such  as  cylindrical  or  annulus)  should  be  solved  in  a  similar  manner  to  the 
present  parallel  plate  problem.  A  potential  source  of  difficulty  in  solving  problems  using 
more  complicated  coordinate  systems  is  the  ability  (or  possible  inability)  of  MAPLE  to 
solve  the  final  equation  for  the  time-averaged  temperature  term. 
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